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ABSTRACT 

The production of ''Li and ^Be during the explosive hydrogen-burning that 
occurs in nova explosions is computed by means of a hydrodynamic code able 
to treat both the accretion and the explosion stages. Large overproduction 
factors with respect to solar abundances are obtained, the exact value depending 
mainly on the chemical composition of the envelope. Although the final ejected 
masses are small, these results indicate that novae can contribute to the ^Li 
enrichment of the interstellar medium. Furthermore, since ^Be decays emitting 
a gamma-ray (478 KeV), with a half-life of 53.3 days, the synthesis of ''Li could 
be tested during the INTEGRAL mission. 

Subject headings: novae, cataclysmic variables — nuclear reactions, 
nucleosynthesis, abundances — gamma rays: theory 
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1. Introduction 

The origin of lithium and other hght elements is still an unsolved problem in 
astrophysics. It is widely accepted that ''Li isotopes are produced during the Big Bang and 
by "spallation" reactions in the interstellar medium by galactic cosmic rays or in flares (see 
[Reeves 1993| for a recent review). Standard Big Bang nucleosynthesis underproduces ^Li 
with respect to solar by more than an order of magnitude (see however the recent paper by 
peliyannis, Boesgaard fc King 19951 ), whereas spallation reactions by galactic cosmic rays, 
produce ''Li and ^Li simultaneously, as well as ^Be, ^°B and ^^B. These two mechanisms are 
unable to account alone for the present ^Li abundance (''Li /H ^ 2x10"^ ). Furthermore, 
they are unable to produce neither the high isotopic ratio ^Li /^Li observed in the solar 
system (''Li /^Li = 12.5 ± 0.2) nor the ^^B/^'^B one (^^B/^^Bsi 4). Recent measurements 
of the lithium isotopic ratio in the interstellar medium ( [Lemoine et al. 1993| , [Meyer et 



al. 1993|) yield values similar to those found in the solar system, indicating that it has 



remained nearly constant or even decreased during the last 4.5-5 Gyr. The contribution 
of a low energy component of the galactic cosmic rays, confined at the source or by stellar 
flares ( [IVleneguzzi, Audouze fc Reeves 1971 , Uanal, Isern fc Sanahuja 1975[ , [Prantzos, (Jasse 



fc Vangioni-l'lam 1993[ ) can account for the boron isotopes but ^Li is still underproduced. 
Therefore, an extra stellar source able to produce this ^Li without generating ^Li has to 
be invoked. The interplay of these sources in the galactic evolution of lithium has been 



extensively studied ([D'Antona fc Matteuci 1991[ , [Abia, Isern fc Canal 1995[) . 



The synthesis of ^Li by a stellar source requires the formation of ^Be which transforms 
into ''Li by an electron capture, being ^Be half-life 53.3 days. As ^Li is very easily destroyed, 
^Be has to be transported to zones cooler than those where it was formed with a time 
scale shorter than its decay time. This beryllium transport mechanism, as first suggested 
by [Cameron 1955[ , requires a dynamic situation, like that encountered in asymptotic giant 
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branch (AGB) stars and novae. Another possibility is the production of hthium and boron 
isotopes by neutrino induced synthesis, during gravitational supernova explosions ( |Woosley 



|et al. 19"90| and [Woosley fc Weaver, 1995]) . The importance of such a mechanism is still a 



matter of debate ( pVIatteuci, D'Antona fc Timmes 1995| ) 



The production of lithium in AGB stars has been extensively studied and these stars 
represent the unique observational evidence of an autogenic stellar origin, since it has been 
observed in them (|Abia et al. 19911 , [Abia, Isern fc Canal 199^ ). The huge abundances of 



lithium found in some AGB stars are a clear proof that these stars are currently injecting 
important quantities of lithium to the interstellar medium. However it is hard to estimate 
their total contribution, since it depends on the estimated number of such stars, that are 
buried by their own wind ( |Abia, Isern fc Canal 199"3| ). 



The production of ''Li in explosive hydrogen burning and, in particular, in accreting 
white dwarfs exploding as classical novae, was first studied with a parametrized one-zone 
model by [Arnould fc N0rgaard 1975| . Later on ^tarrheld et al. 1978| computed the ^Li 



yields by means of a hydrodynamic code. This code simulated the explosive stage of novae 
without considering the accretion phase, i.e., with an initial envelope already in place. 
The conclusion of this work was that, depending on the initial abundance of ^He and on 
the treatment of convection, ''Li could be formed in substantial amounts during explosive 
hydrogen-burning in novae. This problem has been revisited by [Boffin et al. 199"3| . On 
the basis of an extended nuclear reaction network and of updated nuclear reaction rates, 
but adopting again a parametrized one-zone model, they showed that ''Li could only be 
produced in significant amounts at peak densities lower than lO^g.cm"'^, which are lower 
than those predicted by hydrodynamic simulations. They argued that the reason of the 
discrepancy was the neglect of the '^B(p,7)^C reaction in the calculations of [Starrfield et 



|al. 1978| . However, large overproductions of ^Be were found by |Coc et al. 1995| (using the 
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semi-analytical model of |MacDonald 1983| to obtain the temperature and density profiles 



and a complete reaction network that included *B(p,7)^C) showing that the origin of the 
different results was not that reaction. In fact, [Boffin et al. 199"3| also showed by means 



of a two-zone approximation that the efficiency of mixing by convection is a very critical 
parameter, and they stressed the need of a detailed hydrodynamic model to study ''Li 
production more accurately. 

The purpose of this letter is to compute the synthesis of ^Li in both carbon- oxygen 
(CO) and oxygen-neon-magnesium (ONeMg) novae by means of an implicit hydrodynamic 
code, that includes a full reaction network, able to treat both the hydrostatic accretion 
phase and the explosion stage. An estimation of the contribution of novae to galactic 
enrichment is made, on the basis of the overproductions and ejected masses obtained. The 
importance of the initial chemical composition of the envelope is analyzed. 

The detection of ^Li in novae would confirm our theoretical predicition. Furthermore, 
the detection of gamma-ray emission at 478 KeV, corresponding to the decay of ^Be to ^Li 
(half-life 53.3 days) in the early phases of novae by the future mission INTEGRAL would 
also confirm the thermonuclear runaway model for novae and the nucleosynthesis related to 
it. 



2. Model and results 

A one-dimensional, lagrangian, implicit hydrodynamic code has been developed 
following the techniques described in [Kutter fc Sparks 1972| . The code has been built in 
such a way as to enable the study of both the hydrostatic accretion phase and the fully 
hydrodynamic explosion. Detailed nucleosynthesis is obtained by means of an extended 
reaction network, including 100 nuclei, ranging from to ^°Ca, linked through an up to 
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date network including more than 370 nuclear reactions (see Jose 1996 and Jose et al. 1996 



for the details). Concerning the reactions involved in ''Be synthesis, rates are taken from 
the [Caughlan &: Fowler 1988| compilation, [Wagoner 1969| , Pescouvemont 1989| and [Wiescher 



et al. 1989| . Time-dependent convection is included in the code, since the hypothesis 
that the convection time scale is always shorter than the nuclear time scale, inherent to 
time-independent convection, is not always fulfilled. With this method, partial mixing in 
the convective region is included. 

Complete evolution of the accretion and explosion stages of white dwarfs with masses 
ranging from 1 to 1.25 Mq , accreting at a rate of 2xl0~^° M0.yr~^ , with initial luminosity 
1O^^L0 has been computed. We assume that the infalling material is of solar composition, 
but that some mixing process (diffusion, shear mixing) mixes it with the underlying CO (for 
M=l and 1.15 Mq ) or ONeMg (for M=1.15 and 1.25 Mq ) core. This assumption is based 
on the current prediction that enhanced CNO (or ONeMg) abundances are required in order 
to produce a nova outburst and to explain some observed abundances (see P^ivio 199^ for a 



recent review, and Prialnik fc Kovetz 1995| and Politano et al. 1995| , for recent calculations 



of CO and ONeMg novae, respectively). We want to stress that the problem of the initial 
composition of nova envelopes is rather complicated and that it is far from being understood 
in a self-consistent way. Studies of diffusion during accretion onto CO white dwarfs have 
been carried out ( [Kovetz fc Prialnik 1985| and [iben, Fujimoto fc MacDonald 1992[ ), but it 



is not clear if enough enhancements of heavy elements are obtained in the ejecta. For the 
ONeMg white dwarfs, these studies are still lacking. Therefore, a compromise is to adopt 
some percentage of mixing with core abundances. We have adopted a 50% of mixing by 
mass with core abundances, as was done in the work by Politano et al. 1995 . 



The chain of reactions leading to the formation of ^Be has been extensively discussed 
in Boffin et al. 199^ . During hydrogen burning, the formation of ^Be proceeds through 
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^He(a,7)^Be from the initial ^He as (p,7) reactions cannot bridge the A=5 gap. It is 
destroyed by '^Be{p,-f)^B, followed by either SB(/?+)2^He or SB(p,7)9C(/3+,p)2^He. However, 
at high temperature the photodisintegration of ^B, ^B(7,p)^Be, becomes faster than proton 
capture on ^Be. In these conditions, the effective lifetime of ^Be can become larger than 
the time scale of the outburst ([Boffin et al. 1993|) . For typical densities at the base of 
the envelope at the onset of explosion, this would happen only above Tg ~ 1. Below this 
temperature destruction by ^Be(p,7)^B is efficient. Other destruction mechanisms are 
''Be(Q;, 7)^^C and beta decay to ''Li . Radiative alpha capture on ^Be is always slower 
than proton capture as long as photodisintegration of ^B is not efficient. The half- life of 
^Be (53.3 days) would allow the formation of ^Li long after the outburst when the cooler 
envelope would prevent the rapid destruction of this fragile isotope. Since ^Be is more 
efficiently destroyed than produced below Tg ~1 and since it originates only from initial 
■^He , it can only be formed during the outburst if enough ^He survives the initial phase 
when the hydrodynamic time scale is much longer than in the explosive phase. 

The '^He found in the envelope originates from the accreted material and from the 
reaction ^H(p,e~z^)^H followed immediately by ^II(p,7)'^He, which increases slightly the ^He 
abundance in the initial phase of accretion. The two major modes of ^He destruction 
are through ^He(^He,2p)^He and ^He('^He,7)^Be. ^He('*He,7)'^Be is always slower than 
'^He(^IIe,2p)^He, except at lower ^He abundances. As noted by |Boftin et al. 199^ , the 
latter reaction is responsible for the logarithmic dependence of the ''Be yield with respect 
to the initial ^He abundance above X(^He)0. This means that hypothetical higher than 
solar initial ^He abundances, related to enriched secondary star envelopes, do not alter 
dramatically the final ''Be yields. The results presented in this paper are little affected by 
nuclear uncertainties as the rates of the reactions of ^He destruction are precisely known. 

For the ONeMg model with M=1.15 Mq , hereafter called ONe model since magnesium 
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is almost absent (see Dominguez, Ibrnambe fc Isern 1993 and Ritossa, Garcia-Berro fc 



[ben 1996|) , we show the profiles of ^He and ^Be abundances along the envelope for different 



times starting at the beginning of the accretion phase in figure |T]. The corresponding 
temperatures at the base of the envelope are 10 ''K, 2x10''' K, 3x10'' K, 5x10'' K, lO^K, 
and the maximum temperature (2x10^ K). An additional model, for which a considerable 
expansion has already occurred (Rwd > lO^^cm), is also shown. Our results indicate that 
^He is destroyed down to abundances between 10~^ and 10"'', by mass, at the end of the 
accretion phase. More specifically, these abundances correspond to the phase during which 
temperatures are around 10® K, allowing the photodisintegration of '^B to prevent ^Be 
destruction (see discussion above). Therefore, the final ^Be abundances are similar to the 
^He ones at this critical phase. It is important to notice that these values are much higher 
than those at the burning shell (~ 10^^ by mass, see figure |l|), indicating that one-zone 
models are unable to provide correct yields. The average mass fraction of ^Be in the shells 
that will be ejected and, thus, will contribute to interstellar medium enrichment, is around 
10~^. Since all the ^Li finally produced comes from the decay of ''Be and as temperatures 
of our last model are low enough to prevent ^Li destruction, the final ^Li yield corresponds 
to the addition of ^Be and ^Li mass fractions in the ejected shells. The final ejected mass is 
1.9xl0~^ M0, with a mean abundance of '^Li by mass of 6.0x10"'' . Thus, 1.1x10"^^ M© 
of ''Li would be ejected (see tablel). 

Concerning the CO cases, ^He destruction is less pronounced, at the same critical 
phase, and this allows the synthesis of a larger amount of ''Be . The corresponding profiles 
of ^He and ^Be are shown in figure ^, for the M=1.15 Mq case. The mean mass fraction 
of ''Li is 8.2x10"^ (the ejected mass of '^Li is l.lxlO~^° Mq and the total ejected mass 
is 1.3x10"^ Mq). The main reason for the different nucleosynthesis in both cases is that 
for CO novae the presence of ^^C implies that the fast reaction ^^C(p,7)^'^N(/9"'")^'^C is 
dominant during the late accretion phase (whereas the energy production through this 
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reaction is lower in ONe novae, as they almost lack from ^^C). Consequently, the duration 
of the phase prior to maximum temperature is shorter in CO novae, preventing an efficient 
^He destruction and thus leading to a larger final amount of ^Be synthesized. 

As a summary of our results (see table 1, where additional cases are shown), 
overproductions of ''Li with respect to the solar abundances between 100 and 2000 are 
obtained, depending mainly on the chemical composition of the envelope, which is related 
to that of the underlying core. It is hard to estimate its real contribution to the ''Li 
enrichment in the galaxy, since theoretical models systematically produce ejected masses 
smaller than the observed ones. For instance, our models typically eject ~ 10~^ M© while 
the estimated mass of QU Vul 1984 (which has been invoked as a true "neon nova") is 
around 10-^ M© (|S aizar et al. 1992| ). For an overproduction factor as large as 2000 and a 



total ejected masses of ~ 10 ^ Mq (two orders of magnitude lower than observed for QU 
Vul 1984), a nova event would produce ~ 10~^° of ^Li . If we adopt the galactic nova 
rate of Pella Valle fc Livio 1994] (20 yr"^) and an age of the galaxy of ~ 10 years, novae 
should produce at least 20 M0 of ''Li . This quantity is clearly smaller than the estimated 
present content of ''Li in the galaxy, ~ 150 Mq , but, given the uncertainties in the ejected 
mass per event the contribution of novae to the galactic content of ^Li cannot be ruled out 
yet. A complete analysis of ''Li yields by novae and their inclusion in a model of galactic 
evolution is out of the scope of this paper and will be presented elsewhere. 



3. Discussion and conclusions 

Our results confirm that nova explosions can produce significant amounts of ''Li . 
Overproduction factors as large as 2000 are obtained. Our results are quite different 
from those obtained from one-zone models ( [Boffin et al. 1993|) , as the most important 
contribution to ''Li enrichment comes from the external shells, where this element has been 
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transported by convection from the burning shell. 



Comparison with the results of |Coc et al. 1995| and [Folitano et al. 1995| shows that 



the behavior of ^Be abundances can only be correctly predicted if the evolution of ^He 
during the accretion phase is accurately followed. It is also necessary to stress that the final 
results strongly depend on the chemical composition at the onset of the explosion. If the 
underlying white dwarf is a CO one, the ''Li abundances are about one order of magnitude 
larger than if the white dwarf is an ONe one. 

Since the decay of ^Be to ^Li emits a photon with energy 478 KeV, during a phase 
in which the envelope is very transparent, this transition could be detected by the future 
INTEGRAL mission (with a sensitivity around 6x10^^ at this energy). The flux of the 
''Be decay line is: 

F(counts.s-.cm-) = 2.2 x IQ-^^^^e-^^^ 

For an ejected mass of 10^^ M0 , with an abundance by mass of XfBe)=8xlO-^ , the 
'^Be decay line would be detectable just after the outburst only for a nova closer than 0.5 
Kpc. But for an ejected mass of 10~^ or 10~^ M© , more in accordance with observations, 
the lower limit for the distance would be 1.7 or 5.4 Kpc, respectively. This detection 
would provide a confirmation of the theoretical models of novae and also ensure that ''Li is 
produced in these scenarios, encouraging a deep search of this element in novae. 



Research partially supported by the CICYT (ESP95-0091), by the DGICYT 
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CESCA/CEPBA". 
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Table 1: ^Li yields and ejected masses for some nova models 



Comp. 


Mwd(M0 ) 


M(M0.yr-i ) 


X(^Li ) 




N(^Li/H) 
JV(7Li/H)Q 


Ml(Mo ) 


M?Li(Mo ) 


CO 


1.0 


2x10"^° 


3.1x10" 


-6 


742 


2.3x10"^ 


7.1x10-" 


CO 


1.15 


2x10-^0 


8.2x10" 


-6 


1952 


1.3x10-5 


1.1x10-1° 


ONe 


1.15 


2x10-^0 


6.0x10" 


-7 


143 


1.9x10-5 


1.1x10-11 


ONe 


1.25 


2x10-1° 


6.5x10" 


-7 


155 


1.8x10-5 


1.2x10-11 


ONe 


1.25 


2x10"^ 


7.9x10" 


-7 


187 


8.3x10-^ 


6.7x10-12 
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Fig. 1. — Profiles of ^Be (upper panel) and ^He (lower panel) abundances along the envelope 
for different times from the beginning of accretion up to the ejection of the envelope, for 
a I.I5M0 ONe novae accreting at a rate M=2xl0~^° M0.yr~^ . The succesive models 
correspond to temperatures at the base of the envelope 2x10^ (solid), 3x10^ (dot), 5x10^ 
(short dash), 7x10^ (long dash), 10^ (dot - short dash), 2x10^ K (T^ax) (dot - long 
dash) plus an additional case, for which a considerable expansion has already occurred, 
Rwd > lO^^cm (short dash - long dash). The upward arrow indicates the base of the ejected 
shells. 



Fig. 2.- 



Same as previous figure but for a CO nova. 
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